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Abstract 

Azimuthal distributions of particles produced in nucleus-nucleus collisions are mea- 
sured with respect to an estimated reaction plane which, because of finite multiplicity 
fluctuations, differs in general from the true reaction plane. It follows that the mea- 
sured distributions do not coincide with the true ones. I propose a general method 
of reconstructing the Fourier coefficients of the true azimuthal distributions from the 
measured ones. This analysis suggests that the Fourier coefficients are the best ob- 
servables to characterize azimuthal anisotropies because, unlike other observables such 
as the in-plane anisotropy ratio or the squeeze-out ratio, they can be reconstructed 
accurately. 
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A characteristic aspect of collective behavior in nucleus-nucleus collisions is that the di- 
rections of the outgoing particles are correlated to the orientation of the impact parameter [I]: 
azimuthal distributions measured from the reaction plane (which is the plane containing the 
impact parameter and the beam axis) are not uniform. At energies between 100 MeV and 
2 GeV per nucleon, several experiments have measured azimuthal distributions of charged 
particles 0, identified protons and light nuclei || §], charged and neutral || pions, 
neutrons J7|, and A baryons f9j. Similar results have recently become available from the 
ultrarelativistic heavy ion experiments performed at the Brookhaven AGS JT(J and at the 



CERN SPS[|TT| where proton and charged pion azimuthal distribution have been measured. 
In all these analyses, the azimuthal angle is defined with respect to an estimated reaction 
plane which, because of finite multiplicity fluctuations, differs from the true reaction plane. 
It has been recently emphasized that a determination of nuclear equation of state from az- 
imuthal anisotropies requires a good accuracy Jl2]|. It is therefore important to correct the 



errors which are made in determining the reaction plane. 

We propose a systematic procedure to reconstruct the true azimuthal distribution from 
the measured ones. Normalized azimuthal distributions can be expressed as Fourier series: 



dN 




2^c„ cos n< 

n>l 



where c n = (cos ncf)) , the brackets denoting average values, and we assume that the azimuthal 
distributions are symmetric with respect to the reaction plane (i.e. even in <p), which holds for 
spherical nuclei. The knowledge of all the Fourier coefficients (cosn0) allows to reconstruct 
the full distribution using Eq. ([[]). 

In an actual experiment, the reaction plane is not known exactly. It is reconstructed 
event by event from the reaction products. The reconstructed plane differs in general from 
the true reaction plane by an error A(f>, which varies from one event to the other. Thus, the 
measured azimuthal angle ip is related to the true azimuthal angle (f) by ip = <p — A0 (see 
Fig.l). Averaging over many events, assuming that <ft and A<p are statistically independent 
(this assumption will be discussed below), one obtains the following relation between the 
measured and true Fourier coefficients |T^, [L3|, [T^j : 



(cos nip) = (cosn0)(cosnA0). (2) 

From Eq.(^), we can reconstruct the true distribution once the correction factor (cosnA0) 
is known. 

Before we calculate (cosnA0), let us briefly comment Eq.(^). This equation shows that 
the measured anisotropies are always smaller than the true ones: they are smeared by the 
error A<p. More precisely, if the probability distribution of A0 has a typical width 5, i.e. if 
S denotes the typical error made in determining the reaction plane, (cos nA(f>) will decrease 
with n and become small for n > 1/5. This has two consequences. First, the higher order 
coefficients disappear in the measurement procedure. Indeed, all the distributions measured 
so far are well reproduced by keeping only the first two Fourier coefficients n = 1, 2 in Eq.(jl]) 
@? |§! §]) 0' 03 • However, higher order components might be sizeable in the true distributions. 
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The second consequence is that observables involving higher order Fourier coefficients cannot 
be reconstructed accurately. This is unfortunately the case for widely used observables such 
as the in-plane anisotropy f7|, |15| 

_ (dN/dcfi^o , , 

-<~lm - plane — / 7 at / 7 i\ \y) 
(dN/ 0(0)0=180° 

or the squeeze-out ratio [@, || |. §] 

P (dN/d(p) 0=90° + (dN/d(fr) ^=270° . . 

Squeeze - (^/^^ + ( dN / fy)^ W 

which both involve an infinite number of Fourier coefficients (see Eq.(|l])). Furthermore, 
Fourier coefficients, which are integrated quantities, are also easier to evaluate in theoretical 
models (especially in Monte-Carlo calculations) than observables using the value of the 
distribution at a specific point. 

We now recall how the orientation of the reaction plane is estimated from the reaction 
products: in high energy collisions (E/A > 100 MeV), the projectile (target) fragments are 
deflected away from the target (projectile). Therefore, the vector obtained by summing all 
the transverse momenta of the particles produced in the projectile (target) rapidity region 
is parallel (antiparallel) to the impact parameter. More generally, one constructs a vector 
QUI: 

N 

Q = ^2w k u k (5) 

k=l 

where the sum runs over all the detected particles in the event, is the unit vector parallel 
to the transverse momentum of the particle, and w k is a weight which may depend on the 
type of particle, its rapidity and transverse momentum. The choice of w k is to a large extent 
arbitrary. Danielewicz and Odyniec |16| choose w k = px for y > 0.3, w k = —px if y < —0.3 



and w k = if \y\ < 0.3. Many alternative definitions have been used flTTL [TBI, some of which 



do not require particle identification. This method of determining the reaction plane is 
commonly referred to as the transverse momentum method. In experiments which measure 
the energy deposited in a calorimeter, without counting individual particles, an equation 
similar to Eq.(p|) must be used, with w k replaced by the energy deposited in each segment 



of the calorimeter ]19] . We come back to this case at the end of this paper. 

For an ideal system with infinite multiplicity N, Q lies in the true reaction plane, and 
azimuthal distributions can be measured from Q. In an actual experiment, the multiplicity is 
finite, which has two effects. First, because of statistical fluctuations, there is a deviation A0 
between the true reaction plane and Q. Another effect of the finite multiplicity is that when 
one measures the azimuthal angle of a particle with respect to Q, there is a correlation if the 
particle is included in the sum in Eq.(|3). This can be avoided by constructing a new vector 



Q obtained by summation over the N — 1 remaining particles ||16|| . Then, it is reasonable to 
assume, as we have done, that <fi and A0 are statistically independent. However, this is not 
quite true for small nuclei where the constraint of global momentum conservation creates 



important correlations. A method to subtract these correlations is described in |20| . 
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We now proceed to evaluate (cosnA<p). We first show that the distribution of A0 is a 
universal function of a single real parameter x, which measures the accuracy of the reaction 
plane determination, and scales with N like y/~N. It is normalized in such a way that for 
large x-, the standard angular deviation is (A0 2 ) 1 / 2 = l/(xv / 2)- Then we express (cosnA0) 
as a function of x- Finally, we explain how to extract x from the data. 

We consider a large sample of events having the same magnitude of impact parameter. 
Experimentally, this can be done approximately by selecting events having the same multi- 
plicity, or the same transverse energy, or the same energy in a zero degree calorimeter. The 
number N of particles entering the definition of the vector Q in Eq. (|5]) is usually much larger 
than unity. Then the central limit theorem shows that, for given magnitude and orientation 
of the true impact parameter, the fluctuations of Q around its average value, (Q), are gaus- 
sian. Note that both the magnitude and angle of the vector Q fluctuate (see Fig.2). We 
choose the direction of impact parameter as the x-axis. Then Q = Q(e x cos A0 + e y sin A<f>) 
and (Q) = Qe x , and the two dimensional distribution of Q takes the form 

dN 1 ( |Q — <Q)| 2 ^ 1 ( Q 2 +Q 2 ~2QQ C osA(f) \ 

QdQdAcj) ~ trx 2 6XP ^ a 2 ) ~ ttcx 2 6XP ^ a 2 )' W 

We have assumed that the fluctuations are isotropic. This will be justified later. Note that 
Q scales like N while a scales like y/N . 

Eq.© can be easily integrated over Q [^1J to yield the distribution of A0: 

^ = 1 exp(- X 2 ) {l + zJH [1 + erf(z)] exp(^ 2 )} . (7) 

where z = x cos A0 and erf(x) is the error function. This distribution depends on Q and a 
only through the dimensionless parameter x = Q/cr. The Fourier coefficients are most easily 
calculated by integrating Eq.@ first over A<p and then over Q |)H |: 




/ A \ V ^ —V 2 II 

{cosnA(p) = ~^~X e 

where Ik is the modified Bessel function of order k. The variations of the first coefficients 
with x is displayed in Fig. 3. As expected, cosnA0 decreases with increasing n, and becomes 
vanishingly small for n ^> x- Experiments report values of (cos A(f>) ranging from 0.35 [13 



for light nuclei (A ~ 20) to 0.94 |J for the heaviest ones (A ~ 200), corresponding to values 
of x between 0.4 and 2.2 respectively. Fig.3 shows that the corrections are important in this 
range. 

If x ^> 1) the distribution of A<p, Eq.(|7|), becomes approximately gaussian 

dN - ^exp(- X 2 A0 2 ) (9) 



dA4> v/tt 
and the Fourier coefficients are given by 

(cosnA0) ~ exp(-n 2 /4x 2 ). (10) 
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In the limit % -C 1, the n th Fourier coefficient is of order \ r ' 



(cos nAcfi) 



71 



2 n T (2±i 



where T is the Euler function. 



We now turn to the determination of x- The most widely used method |J 0, [H| to 
estimate the accuracy of the reaction plane determination is to divide each event randomly 
into two subevents containing half of the particles each, and to construct Q for the two 
subevents [|1(J. One thus obtains two vectors Qj and Q H . The distributions of Qj and Qjj 
are given by an equation similar to Eq.(|6]). However, since each subevent contains only N/2 
particles, the corresponding average value and fluctuations must be scaled: Qi = Qn = 0/2, 
°7 = &ii = 0"/v2j an d therefore xi — Xn — Xl The distribution of the relative angle 



A(f)R = \A(f)j — A(j)jj\ can be calculated analytically (see the Appendix of and the note 
added in proof) 



diV 



dA0 



R 



7T 



1 + X 2 i) + z [I (z) + LoC*)] + x) lh(z) + U{z)] 



(12) 



where z — Xj cos ^■'Pr an d L and L x are modified Struve functions |22j . This distribution 
is normalized to unity between and tt. 



The value of x can be obtained by fitting Eq.(|i~2"D to the measured distribution. However, 
it can be calculated more simply from the fraction of events for which > 90°, which is 
calculated by integrating Eq.(|i~2"D over AcpR. 



iV(90° < A<p R < 180°) _ exp(-Xj) _ exp(-x 2 /2) 
A^(0° < A<p R < 180°) ~ 2 ~ 2 ' 



(13) 



Alternatively, one can obtain x by measuring [19, 21 



(cosA(j) R ) = (cos A0/)(cos A<j>u) 



* 2 - 
aX e 



-X 2 /2 



J (x 2 /4) + h (x 2 /4 



(14) 



where we have used Eq.(§) with n = 1 and x replaced by xi — x/V%- The variation of 
(cosA0a) with x is displayed in Fig.4. Eq . (|14]) is more reliable than Eq.(^) if x is large, 
for in this case the ratio in Eq.(|T^) is very small and is therefore subject to relatively large 
statistical fluctuations. Other methods to measure x are described in ||2l| , p3| . 

We finally justify the hypothesis that was made in writing Eq.(||), namely that the 
fluctuations have the same magnitude in both x and y directions. This is not a consequence 
of the central limit theorem, which only ensures that the two dimensional distribution of Q 
is gaussian. The most general gaussian compatible with the symmetry A0 — > — A0 can be 
written as 



dN 
QdQdA(f) 



1 



na T a. 



exp 



f) 



(QcosA0-Q) 2 Q 2 sin 2 A0 



<Jt, 



(15) 
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The quantities a x and a y characterize the fluctuations along the x and y axis, which are 
a priori different (the circle in Fig. 2 should then be replaced by an ellipse). Q, a x and o~ y 
are related to the azimuthal distribution of particles in the following way. Assuming for 
simplicity that the multiplicity N is the same for all events in the sample (it is at least 
approximately true since the impact parameter is fixed), we get from Eq.fl5|) 



Q = (Q • e^) = N(w cos i 



(16) 



where <ft is the true azimuthal angle of particles, and the last average involves all the detected 
particles of all events. Similarly, the fluctuations in the x and y directions are given by 



((Q-e z ) 2 )-g 2 J = 2N 
ct 2 = 2((Q-e y ) 2 ) = 2N(w 2 sm 2 > 

We define the average fluctuation a by 
o 1 



(w 2 cos 2 0) — (w COS ( 



(17) 



a' 



{<?) - Q 



N 



(w 2 ) — (wCOS< 



and the anisotropy of the fluctuations by 

1 



{pi - °i) 



N 



(w 2 cos 2(f)) — (wcos< 



(19) 



Three cases must be distinguished, depending on the relative magnitudes of Q, a x and o y . 
• (A) Azimuthal anisotropies are small, i.e. the Fourier coefficients of the azimuthal 



distribution are much smaller than unity. Then (w cos 2(f)) and 



[WCOSt 



are both 



18Q and (pj|) that a x ~ o y ~ a: our 



small compared to (w 2 ). We deduce from Eqs. 
assumption is justified in this case. 

(B) There are situations where (w cos 0) 2 cannot be neglected compared to (w 2 ). This is 
the situation when the flow is strong. Using Eqs.fll6|) and ( |PTD and the fact that N ^> 1, 
we see that in this case Q 3> o~ x , a y \ fluctuations are small and A0 ~ Q-e y /Q 1 (see 
Fig. 2). Therefore, the distribution of A0 involves only a y , not a x . Although a x and 
<j y may differ in this case, one can replace a x by a y without altering the distribution 
of A0. Note that in this case, x = (Q/a ^> 1, hence the distribution of A<p reduces to 
its asymptotic form, Eq.(|). 

(C) Finally, there is the case when (w 2 cos 20) is of order (w 2 ) while (it; cos 0) 2 is much 
smaller. In this case, there is a strong anisotropy in the second Fourier component, 
which should be used to determine the reaction plane. Instead of Q, one should 
construct the vector Q 2 defined as 



N 



Q 2 



fc=i 



ifc + e y sin 2(j> k ) 



(20) 



with the same notations as in Eq.(^[), and w' k is an appropriate weight. The azimuthal 
angle of the reaction plane is estimated as half the azimuthal angle of Q 2 , hence it is 
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defined modulo 7r, i.e. one cannot distinguish </> and + n. This method is equivalent 
to the diagonalisation of the transverse sphericity tensor, which has been claimed 
to be more efficient than the transverse momentum method at intermediate energies 
(E/A < 100 MeV, see @ and in particular Eq.(12) of P5fl ) and at ultrarelativistic 



energies where the flow angle is very small due to increasing nuclear transparency 



21], pq| . Then, the procedure described in this paper can be applied to reconstruct 
the azimuthal distributions, replacing Q by Q2. However, since the azimuthal angle is 
defined modulo tt, only the even Fourier components can be reconstructed. 

Before we conclude, let us briefly comment on corrections which have been applied, in 
previous works, to the measured azimuthal distributions. The method proposed by Tsang 
et al. |25|, [27J is based on an ansatz for the distributions of and A<f>, which are assumed 
to be proportional to exp(— u 2 sin 2 (ft), where u is a fitted parameter. Our analysis is more 
general in the sense that it does not make any a priori hypothesis on the shape of azimuthal 
distributions. The analysis of Demoulins et al. |D| is similar to ours, although limited 
to n — 1 and n = 2. The main difference is that they calculate (cosA0) and (cos2A</>) 
independently from measured quantities. The analysis done by the E877 collaboration |TU| 
is similar, in the sense that the corrections (cosnA0) are extracted directly, for each value 
of n, from measured correlations between subevents. On the contrary, our method allows to 
calculate simply all the Fourier coefficients (cosnA0) from Eq.flHD or Fig.3, as soon as we 
know the first one, (cosA0). 

In conclusion, we have described a simple procedure to reconstruct the true azimuthal 
distributions from the measured ones by means of analytical formulae. Let us summarize this 
procedure: given a large sample of events in a restricted centrality interval, one measures the 
distribution of the relative angle between subevents A0ij, from which one extracts, using 
Eq.([I3]) or Eq. fl!H|) , the crucial parameter \ which measures the accuracy of the reaction 
plane determination. Then one uses Eqs.@ and @ to reconstruct the Fourier coefficients 
of the true azimuthal distribution from the measured ones. 

A first assumption, on which our calculation relies crucially, is that the fluctuations of 
the momentum transfer Q are gaussian. This is reasonable only if the multiplicity is large 
enough. The method should not be applied if, for instance, a single big fragment of the 
projectile gives the largest contribution to Q. In order to test the validity of the gaussian 
hypothesis, one may check that Eq.(|T2"D provides a good fit to the measured distribution of 
the angle A0#. Alternative tests are proposed in [23]. Our second assumption was that 



the subevents used in determining the value of x are equivalent. This holds strictly only 
if the subevents are defined by a random selection of N/2 particles. The subevents are 
not equivalent if they correspond to different pseudorapidity windows as in [TIJ. However, 



Eq.(ll) can still be used to calculate the correction for arbitrary n, once the correction for 
n — 1 is known. 

We have introduced a parameter \ to measure the accuracy of the reaction plane deter- 
mination. Unlike other quantities frequently used in the literature such as (cosA0) ||, |13|1 , 
or (A0 2 ) 1 / 2 [|, H [7], [|, X scales simply with the multiplicity like yN . It has a simple physical 
interpretation, being the ratio of the average value of the flow vector Q to the typical statis- 
tical fluctuation a. And furthermore, it can be directly deduced from measured quantities 
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by simple expressions such as Eq.([T 

Our analysis suggests that the Fourier coefficients of the azimuthal distribution, which 
can be reconstructed accurately and are also easy to estimate in theoretical models, are 
the best observables to characterize azimuthal anisotropies. Note that we have chosen to 
measure azimuthal angles around the beam axis. It was argued in Ref.p8| that azimuthal 



distributions should rather be measured around the flow axis, determined by a sphericity 
tensor analysis |29j . But the fluctuations of the sphericity tensor, which is a 3 x 3 matrix, 
are much more complex []30| than those of Q, and cannot be described in terms of a single 
parameter x- There exists no simple procedure to subtract statistical errors in this case. 

What accuracy can be attained with our method? If the azimuthal angles of M particles 
are measured (summing over all events), the corresponding statistical error on the measured 
Fourier coefficient is 1/ v2AT. Systematic errors can be removed, at least partly, using 
a mixed event technique. The error on the true Fourier coefficient will be larger by a 
factor l/(cosnA0), according to Eq.(|2]). However, with high enough statistics, accurate 
measurements are possible even if the reaction plane is poorly known, as is the case with 
light projectiles |13|]. In experiments using heavy ions, it would be interesting to try to 
measure higher order Fourier coefficients with n > 3 which, although probably small, could 
be measured accurately. 
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Figure captions 



Fig. 1: Schematic picture of a semi-central nucleus-nucleus collision viewed in the transverse 
plane (the beam axis is orthogonal to the figure), b is the impact parameter oriented from 
the target to the projectile and Q is the vector defined by Eq.(Q). A particle is emitted along 
the dashed arrow. Its azimuthal angle measured with respect to Q is if), while the "true" 
azimuthal angle is ip + A0. 

Fig. 2: Schematic picture of the distribution of Q, given by Eq.(|6|). The solid thick arrow 
indicates the average value (Q) = Qe x , which lies along the direction of the true impact 
parameter. Q fluctuates around this average value with a standard deviation a, so that a 
typical value of Q (dotted arrow) lies within the dotted circle with radius a. It is obvious 
from this figure that the typical magnitude of A0 is a/Q = 1/x- 

Fig. 3: Solid lines: variation of (cosnA0) with the parameter x, calculated from Eq.([8]). 
The curves are labeled by the value of n. The dotted curves and dash-dotted curves are the 
asymptotic forms given respectively by Eqs.([TT|) and (flO|) . 

Fig. 4: Variation of (cos A0^) with x, given by Eq.flH|). 
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